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Abstract 

A path-integral method effective beyond the perturbation expansion approach is 
suggested to consider the quartic anharmonicity in different spatial dimensions. 
Due to an optimal representation of the partition function, the leading term 
has already taken into account the correct strong-coupling behaviour. In the 
simplest cases of zero and one dimension we have obtained reasonable results 
in a simple way. Then, this technique is applied to the superrenormalizable 
scalar theory <$>\ in two dimensions. This results in an accurate estimation of the 
ground-state energy that provides exact weak- and strong-coupling behaviour 
already in the leading-order approximation. The next-to-leading terms give rise 
in insignificant corrections. 

I. INTRODUCTION 

The anharmonic oscillator plays an important role by providing a theoretical laboratory for 
the examination of new calculational schemes and approximation techniques (see, e.g. 0,0]) in 
quantum mechanics. Besides, it describes more or less real self interaction disturbing the ide- 
alized picture of the free harmonic oscillators. The most known version is the quartic oscillator 
[]3],^[]. On the other hand, scalar fields play a fundamental role in the unified theories of strong, 
electromagnetic and weak interactions. The mechanism of symmetry rearrangement specific 
for scalar fields occurs in realistic scalar field theories [H. At the same time, the mathematical 
structure of scalar field theory is simpler than that of the vector and spinor cases. A native 
and well-known generalization of the quartic oscillator is the model of scalar self-interaction 
<f) 4 in quantum field theory. Particularly, the quantum mechanical quartic oscillator with a 
Hamiltonian 

ff = i(p? + my + 0V) (1) 

can be considered || a theory of a scalar field in one-dimensional (q(t), t e R 1 ) Euclidean 
space-time described by the Lagrangian: 
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L[qA] = \[q 2 (t) 



(2) 



Thus, the scalar field model (0) represents an attractive, simple but nontrivial object for re- 
search. The problem of finding the energy spectrum {E n } and the proper wave functions {^ n } 
of the Hamiltonian can be solved by using the partition function and path-integral formal- 
ism. Note, the corresponding Schrodinger equation cannot be solved explicitly for general case, 
except the weak- and strong-coupling limits. 

Further extension of the model fl2|) can be done by considering the quartic self interaction 
4>\{x) in higher spatial dimensions x G R d , d > 1. Earlier investigations of the triviality 
problem have shown that for dimensions d > 4 the scalar theory of self interaction behaved 
either unstable or trivial. Much interest is attracted to the cases d = 2 and d = 3 demonstrating 
nontrivial phase restructure. 

Nowadays, the anharmonic oscillator is a system well understood in both perturbative and 
non-perturbative aspects |J[7||. 

Typically, the problem is treated by using various perturbative expansion scheme. Par- 
ticularly, the conventional loop expansion [Q gives asymptotic Rayleigh- Schrodinger series for 
energy levels which numerically converges only for j < 1. The loop expansion for composite 
operators || and optimized expansion pT0| , [TT|] provide better convergence for numerical series; 
however, in this case the approximation becomes worse as g increases. To obtain a reasonable 
result, one needs in various resummation techniques for these divergent series [|EJ]. Particularly, 
the Heaviside transformation of the mass parameter [13 and a modified Laplace transformation 
method [TJJ can result in good approximations to the ground state energy, but the perturbation 
series depends on cut-off parameter x* and converges very slowly [14]. The Pade approximants 



and the Borel transform can be used to summarize the formally diverging series But, rig- 



orous proofs of the convergence can be made only in particular cases and the convergence 
is nonacceptably slow. 

Going beyond the perturbation formalism, a variational interpolational method W7\ and 
the strong-coupling expansion [jrj| are used to estimate the energy spectrum of the AHO. 
Also, by applying the renorm-group equation, one obtains a resummed perturbation series [|TTfl. 
The renorm-group methods may improve the perturbative expansions fpOJpJJ. However 



the 

agreement with the WKB result becomes worse in the higher orders than the fourth at which 



the agreement is the best. Recently, it is have shown that multiple-scale [22] and reductive 
perturbation theory [12| can be successfully applied to the quantum AHO to construct the 
asymptotic behaviour of the wave function for large x. A similar method had been proposed 
in ||23|| . But, it is not trivial to identify the secular terms for the wave functions which can be 
made to vanish at a renormalization point for general case. 



Within the method of orthogonal polynomials it has been established [25] that the solutions 
of the bound state equation for the quartic anharmonic oscillator are related to a sequence 
where the potential is a even polynomial of degree six (a quasi-exactly solvable model). A 
regular iterative method for the ground state energy of quartic oscillator as an analytic func- 
tion of coupling constant [2JJ is formulated. But this method introduces a PI averaging of the 
complicated many-point potential of the logarithmic type that makes it hard to obtain explicit 
solutions. By using an extension of quasi-exactly solvable models [26| to Bose systems, one 
obtains (in some particular cases) simple exact expressions for several energy levels of an an- 



harmonic Bose oscillator |27[] beyond the exploit perturbation theory. A relationship between 
the x 2M AHOs and the A 2 m-i TBA systems have been conjectured [28| by using an alternative 
integral expression for the spectral determinants of the quartic oscillator. 
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Therefore, it is important to build a method which is effective in the strong-coupling. On 
the other hand, it has to work even in the case of renormalizable theories (e.g., (p^i 00 )) with 
formally divergent or complex functionals, unlike the variational methods. Besides, the desired 
method should be relatively plain and has to catch the main strong-coupling contribution within 
the leading-order approximation. Below we represent a path-integral technique obeying these 
features. 

In our earlier works [p9"[|3~0f the problem of strong coupling regime and phase structure 
of the y>\ quantum field theory has been investigated. We have found that an approximate 
solution of the problem can be obtained by combining the canonical transformations of field 
variables and normal ordering of the creation and annihilation operators. In the lowest ap- 
proximation, the solution obtained there is identical to the result of the variational method 
of the Gaussian effective potential [[31]. Meanwhile, in contrast with the variational method 



our suggested technique ensures a canonical structure of the theory and, hence, provides cal- 
culation of perturbation corrections over effective coupling constants which can also be used to 
control an accuracy of approximation. Later on, the method was generalized and applied to 
renormalizable quantum field theories, path integrals and quantum mechanics of bound states 



32] . 



In the present paper we demonstrate a significant modification of our method suitable to 
investigate the 4 model effectively even for large coupling. Below we consider the case of low 
dimensions, namely, d = 0,1,2 to avoid nonprincipal but lengthy renormalization procedures 
occurring for d = 4 — e, e — > 0. In fact, an extension of our method from d = 2 to d = 4 — e does 
not meet any principal troubles, because only renormalization of the mass is sufficient to remove 
all divergencies in both models. Therefore, upgrading the dimensionality to d = 4 — e will affect 
only some coefficients and factors, not changing the very form of our final expressions. However, 
we should note that the consideration of the theory <p\ requires a serious modification of our 
method because one has to workout first an effective technique to perform the renormalization 
of the coupling constant. 

Our basic idea is to find a optimal representation of the partition function (e.g., see (|i~7D) in 
which the main contribution is taken into account correctly for large g 2 already in the lowest 
approximation. For this purpose we will use the Gaussian representation of the initial quartic 
interaction in the form of an integrable quadratic one by involving an auxiliary integration as 
follows: 



, e ^l~{qD- l q)-^-q 4 \ oc JdAexp l~ - ^trln(l + 2igqD )\ . (3) 



Then, the inner path integral over can be taken out explicitly and the remaining integral 
over auxiliary variable (field) A can be evaluated effectively, in particular, by isolating the total 
Gaussian contribution. 



II. A NON-GAUSSIAN INTEGRAL: D = 

Below we demonstrate the basic idea of our approach on the simplest case of a plain non- 
Gaussian integral. 

First, we introduce a notation indicating an averaging with respect to a normalized Gaussian 
measure as follows: 
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<(•)>, 



dx_ 2/2 



2tt 



0, (i> x = i. 



(4) 



Let us define a "partition function": 



J] = ( exp <j - yx 4 + ixJ 



(5) 



Then, we can consider the zero-dimensional analogues of the Quantum Mechanical ground- 
state energy and Green function as follows: 



exp{~E(g)} = Z[g,0], G{g) 



1 d 2 Z[g,J] 



Z[g,0] dJ 2 



(6) 



j=o 



As is known, these integrals are not solvable exactly. According to our approach we write 



Z[g, J] = (^ J dxexp j — ^ — igx 2 A + irrJ j 



=( exp 



In (1 + 2igA) + 



J 2 



1 + 2igA 



(7) 



To build the optimal oscillator basis, we have to eliminate the linear A terms out of the 
exponent as follows: 



aA iqA 
— = 



g —a(l + a), h = — = 



g 2 a(l + a) 



2ig 1 • a u 2 

After a shift A — > A + a/2ig with a parameter a > we obtain 



ig Ag 2 v 7 l + a + 2igA 



Z[g,J] = (^expj- 
First we rewrite the logarithm in the exponent as 



2iqA 2q 2 A 2 

In (1 + a + 2igA) = In (1 + a) + + — ^ — + ln 2 (1 + a + 2igA) , 

1 + a (1 + a) 



where 



(8) 



(9) 



(10) 



ln 2 (l + z) = In (1 + z) -z + 



VRez > -1 



Define a function 



W(s) = --m 2 (l + isA) , A = 



2a 



1 + 2a 



Then, we obtain 



Z[g, J] = exp <J — 



a 2 ln(l + a) 
8? 



exp < W(s) 



< I. 



J 2 



2(1 + a) (I + is\) 



(12) 



(13) 



Accordingly, 

Systematically expanding the exponentials in ( fH] ) we estimate the "ground-state energy" 

N 

E N (g) = E (g) + £ AE^g) , AE x (g) = - (W(s)) s (15) 



and the "Green function" 



1/1 



~[ l + a\ l + «sA/s 

where E (g) and GoCfl 1 ) are the leading-order approximations. 

In Table 1 we represent our estimates due to (15) and ([16|) compared with the exact nu- 
merical solutions to (H). We see that our approach describes well these functions even within 
the zero-order approximation. The next-to-leading corrections systematically improve the ob- 
tained results approaching close to the exact solutions. Our technique works well at least in 
zero dimension. 



III. ANHARMONIC OSCILLATOR: D = 1 

The most known example of the quartic interaction is the anharmonic oscillator g 2 q 4 /2 in 
Quantum Mechanics. The main problem is to calculate its energy spectrum and to find the 
corresponding wave functions. 

Conventionally, the problem is treated by using various perturbative methods (e.g., ||[|[IT|). 
However, the ground-state energy in the naive perturbation series || reads, E = m J2%Lo A n • 
(g 2 /m 3 ) n , where the coefficient A n grows as A n ~ (— 3/2) n r(rz) for large n ||33|| . Thus, the 
series formally diverges for any finite coupling. Only for sufficiently small coupling constant, 
g 2 /m 3 ~ 0.1, the series becomes numerically useful by the appropriate truncation. Beyond 
the perturbative approach, variational [[IT]] and strong-coupling expansion methods Jl8| can be 
applied to this problem. 

Our effective nonvariational approach described in previous section can be extended to the 
quantum mechanical anharmonic oscillator as follows. 

Consider a normalized partition function 

Z v (g) = J Dq e-K^o ^-^fe*) , Z v (0) = 1 , 

where q(t) is the position vector and t G [— T, T], Jdt = V — » oo. 

The differential operator Dq 1 and its Green function D obeying the conventional boundary 
conditions are: 

s) = (~ + l) 5(t - s), D (t) = D{k) = . 
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For further convenience we use the unit mass m — 1. 

First, we introduce an additional integration variable (an auxiliary trajectory) <fi to convert 
the initial quartic interaction to an integrable quadratic one by using the following Feynman 
representation 

Z v {g) = J D(j> e-^W j Dq e -\^-\)- ig { q ^) = j D<f> e -iW)-^rin[i+2^A,] _ (17) 

After this transformation, the interactional functional reaches its minimum distinct from 
the origin = 0. Therefore, a shift of the variable may be suggested 



a 



Then, one rewrite 



Z v (g) = exp \ — (aa) - -Trln[l + aD 



— iaa) ' 

8g 2K J 2 

fDcf> exp |-^(#) - ^-(a0) - ^Trln[l + 2ig<f>D}\ , 



2ig 
where 

1 rdk e ikt e^l'l , 

D = j— — FT A) > D(t )= 7TTT- -l = -Z—> /i = v / TT^. 
1 + aV J 2tt k -\- /j, z 2/a 

The optimal value of the shift parameter a obeys the equation: 

4 

In particularly, a = g$ for g — > oo. By using the optimal shift we rewrite 



Zy(sO = ex P \ -^( aa ) ~ o Tr M 1 + aD o] 



J exp|-i(0P-V)-^Trln 2 [l + 2i# J D]} , 



where 



= 5(f) + 2g 2 D 2 (t) , P- X (A: 2 ) = 1 + 



k 2 + 4(1 + a) ' 



p(t) = 5(f) _ ^ e — 1*1 , V(k 2 ) = 1 - , v = v^Te^ . 

v k z + v z 

By introducing a normalized Gaussian measure 



VdeW 

we obtain 
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Z v (g) = e -V(Eo( 9 )+AE( g) ) ; 

E (g) = i {-^ + ~Trln[l + oD ] + ~Trm[l + 2^/J 2 ] 



1 a 1 / 3 

-V 1 + a + W 1 H — a . 



2 Sv/TT^ 2 

A£(s) = d£ c exp{~Trln 2 [l + 2i# J D]} . (18) 

We consider E (g) the leading-order approximation to -E'(g). It can be compared (in Table 
2) with our earlier result obtained within the oscillator approximation method 



E (a) = -1 + m U r(2 - ^ + + 3P) + a ^~ aT ^ ± 5 ^ 1 (19) 

E osc {g) 2 + Tp \8pr(i + p) + 6er(i + P ) + ioe r(i + P )J 1 j 

and with a direct numerical calculation E num (g) (see, e.g. 0). 
In the weak-coupling limit (g <C 1) we have: 

%)=^ 2 - + o(s 6 ) , *Us) = \a 2 - ^X™/ + °^ ' 
^)=^ 2 -|^ 4 + o(/). (20) 

For g —>■ oo one obtains 

Eoig)^ 1 \Jl~fj= 9 1 0.599745 , E osc {g) - <?§ 0.531248 , 
E{g)^g% 0.530181. 

We observe that our leading term E (g) provides result about 15 percent worse than E osc (g), 
but our present technique with only parameter a is much easier than the oscillator approxima- 
tion with two variational parameters and complicated transformations. 



A. Next-to-leading Correction 

Unfortunately, we cannot calculate explicitly the last path integral in (0). On the other 
hand, we have already factorized out the main (generalized Gaussian) contribution e~ VE °^ 
and therefore, the remaining part e~ VAE ^ should not result in a relatively strong correction. 
Therefore, we develop a systematic scheme to estimate the non-Gaussian correction. 

We define the next-to-leading non-Gaussian term as follows: 

A£i( 5 ) = ^JdZ v Tr (ln[l + 2ig<f>D] - 2g 2 [<J ) D ( f ) D}) . 
Taking into account the symmetry <-> — we rewrite 

AE^g) = -L JdZ v Tr (ln[l + Ag^D^D] - ig^D^D]) . 
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Then, going to a new scale 

t-*-=t, k^^/ak, b(k)^—— — -, V(k) -> 1 



y/a ' a(fc 2 + l + l/a)' v/ A; 2 + 6 + 4/a 

we rewrite 

2 

AE 1 (g) = A- Tr | dE c (ln[l + 0] - 6) , = . 

In general, any direct evaluation of AE^g) represents a heavy and lengthy procedure. But 
we can easily estimate it as follows. First, we note that function ln(l + x) — x is concave while 
ln(l + x) — x + x 2 /2 is a convex one. Then, by using these properties, we can easily find a lower 
and an upper bound to AE\(g) as follows: 

AE^(g) < AE]_{g) < AEf(g) , 

^)4W + «]-(e)-M], 



AVa 

,2 



AjBl+( ^ = ik {ln[1 + (0)] " (0)} ' (( - }) = V I d ^M-) • 

We calculate 

^ = 4 1 2^ ( l ~ p 2 + 6 + 4/a) I to k 2 + l + ifa (k + p) 2 + 1 + 1/a = 4Q(a) ' 

„, , 2V1 + a + V4 + 6a 

Q(a) = a 



4(1 + a)V2(2 + 3a + ^(1 + a)(4 + 6a)) ' 
2 ^ = 16 [F(a) + 2i2(a)] , 

3/2 v/TTa (24 + 31a) - V4 + 6a (12 + 11a) 



16(1 + a)(2 + 3a)(3 + 5a) 



F(a) = a 



V(l + a)(4 + 6a) -4- 6a 
W 32(l + a) 3 / 2 (2 + 3a)(3 + 5a) ' 

For strong-coupling regime g > 1 we obtain 



so 



and 



while 



(0) = jl = 0.816497 , (0 2 \ = - - ^ = 1.62407 
w V 3 \ / 3 115 



-# 2/3 0.139072 < A^i(^) < -# 2/3 0.054897 



# 2/3 0.460673 < E^g) < g 2/3 0.544848 



E (g) = g 2/3 • 0.599745 , E(g) = g 2/3 • 0.530181 . 



We see that the next-to- leading correction AEi(g) is negative and it may lower the energy 
E (g) about 10 percent approaching close to the exact numerical result. 
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IV. SCALAR FIELD MODEL: D = 2 



The scalar ip 4 theory in two spatial dimensions has been intensively investigated f34| [36 



as a simple, but nontrivial example, within which the vacuum exhibits a nontrivial structure 



Within the framework of constructive QFT |40|j3^j4l|j35|j4^] a set of general theorems has 



been proven to establish the existence of nontrivial two-dimensional theories of self-coupling 
scalar field. Unfortunately constructive quantum field theory gave no effective instrument (like 
Feynman diagrams) for the calculation of important physical characteristics of the QFT models. 

An attractive approach to the problem under discussion is the variational method of the 
Gaussian effective potential [^-|46[|. Original investigations in the same direction were made 
in |p7|-[50||. An attempt to go beyond the Gaussian approximation has been made in |51 . 

However, specific features of the variational approach in QFT make their results unreliable 
if the theory has divergencies (for d > 1) in the highest perturbation orders (see, e.g. ||52|| ). 

In this paper we calculate the ground-state (vacuum) energy. We demonstrate an effective 
scheme resulting in the explicit and exact asymptotics of the self energy for strong coupling. 

Consider a superrenormalizable scalar model 02 m two-dimensional Euclidean space-time 
x G V C R 2 . The Lagrangian is given: 

£ = i(0[D-m 2 ]0)-^:0 4 :, (21) 

The kernel □ — m 2 and its Green function are 
D~\x - x') = (-□ + m 2 )5(x - x'), 

^/ s r dk e ikx 7 d(3 _n m 2_£. _ . , 7 dB b 

The super-renormalizable theory 02 contains divergences which can be removed by the 
renormalization of mass and vacuum energy. However both these divergences can be removed 
if the interaction Lagrangian in fl2T|) is written in the normal form: 

: 4 : = 4 - QD (f) 2 + 3D 2 , = (0 2 - 3D ) 2 - 6£> 2 . (22) 

This form of interaction removes all the divergences in this theory. For intermediate calculations 
we shall use the dimensional regularization, i.e. we consider the theory in the space R d (d < 2) 

r d d k e ikx 1 7 d/3 -i m 2 -^ 
D ^ X) ~ J (2^p k 2 + m 2 ~2J (2irf3) d / 2 6 * ' ' 
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2 J (2vr/5) rf / 2 2(27^/2 \m 



o 



All divergences are explicitly removed in final formulae and then, we put the true spatial 
dimension d = 2. For simplicity we omit the index "reg" throughout the text having in mind 
that the regularization has been introduced. 

It is convenient to rewrite the Lagrangian as follows: 
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C = \{<f>[U - m 2 ]0) - t- (f - 3D ) 2 + 3g 2 Dl . 



(23) 



In this paper we shall calculate the vacuum energy. Let us consider the partition function 
and use the Gaussian representation: 



Z v [g] = det(-D +m 



J D$ det(-D + m 2 ) J 5(j) e -^[-°WW-^[^-3D ])+3 5 2 D^ 

J D$ e ^9D (<f)-^ r Trln(l+i2 9 'fD)+39 2 D 2 V ^ 

J D<S> e i 9Do($)-±Trhi(l~4g 2 D D+i2g<i>D)- 9 2 D 2 V ? = e ~ 5 



with normalization Zy[0] = 1. 

Let us consider the interaction part 

U[$\= Trln [l - 4g 2 D D + ilg^D 



Trln 



D- 1 - 4g 2 D + i2g$ 



D- 



e 2 



f (-D+m 2 ) _ _s(_ a+m 2_4 s 2 Do+i2 9 *(x)) 



5 d (x - x') 



(24) 



= J dx[W m 2(x,x\0) - W m 2_4g2 Do (x,x\$)] . 
where the function W is defined as: 

W^fox'IS) = / — e -f (-°-Hn 2 -H2fl*(x)) . _ ^ 

J a 
o 

It satisfies the condition 

W m 2 (x,x'\$ + igA) = W m 2 ^ 2g 2 A (x,x'\<5>). 
The function W can be represented in the form of a functional integral (see Appendix A): 



W m 2(x,x'\$) = 
where G K d and 



da 



-e 2 



_ ifl j< & ( x | H v(i-§)+£(/?)) 
2° / da[Q e 



(25) 



-i f 5 2 (/3) 

with the boundary conditions £(0) = £(«) = 0. 
The Green function of the kernel 



/ da[i] = 1 . 



dp 



; 5((3 - (3') 
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reads 

K([3,[3>) = -h{3-[3'\ + kP + P')- — 
2 2 a 

and satisfies the boundary conditions 

K(P, p') = K((3', 0), K(P, 0) = K(0, a) = 0. 

Then, 

oo 

W m ,(x,x\$) = jDa J da[£] e"^) , J Da . (...) = J e"^ • (...) . (26) 

where 

a 

Bz(y) =Jdp5(y-x- £(0)), J dy B 6 (y) = a . 
o 

so that 

a 

jd(3$(x + m) = Jdy *(y)Bz(y) = ($^) , Jdx(^B ( ) = «($) . 
o 

Finally we get 

U[$] = Jdx JDa Jda[£] [l - e WA,-* 9 (*B s )" _ 

Following the idea mentioned in the previous sections, we go to a shift 

<&(x) — > <&(z) - igA 

where A is a constant depending on g. 
We rewrite 

Zy[g] = e 9 2 [hA 2 +D„A-D*]v J m e ig{D 0+ A)^)-\U^-igA\ 



and 



oo 



U[*- igA] =fdxf 2 e~^~ ■ f da® ■ [l - e^W)-*^)] (2?) 
J J (27r) d / 2 a 1+ 2 J 1 J 

Now we introduce the normal-ordered form: 

e -ifl(*B 4 ) =; e -ip(*B € ) :$ e _£(B 4 B 4 ) ^ 



a 



obeying the following relations 

j m , e -i 9 {*B,) :$= x | d(T[ ^ . e ik(m)-m)) :?= i . 

One gets 
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a 

(B 6 B e ) = J dyB e {y)B € {y) = JJ dfcdfc 8 d {^x) ~ £(&)) 



o 

ddk : e iHm)-m)) . e -^-F a Qih-02\) 



(2vr) c 



((BsBs)) £ +:W a [$ : ? 



where 



<(B € S e )> € = / d*[€\(B&) = JJ d/hdfh J ^e-^^IA-ftl) 

d/% - /3) 



(2vr) d / 2 



and the interaction functional is introduced 



/3 / (2vr) d / 2 



d d 
2' ~ 2 







d d k 



ik(HWi)-aM) . p -¥-F a {\p x -fc\) 



k z s 



{2-kY ' 



It is easy to check that 



e g 2 a{2D -A)-ig{<$>Bz) = e -g 2 aA+2g' 2 aD - S- (B f B 6 ) . . e ~ig(0B e ) 



e -g 2 aA-g 2 N(a) . g-%-:W«K]: 4 



: e^*^ :„ +1 - ig{<S>B,) 



where 



-g 2 N(a)=2ag 2 D -^((BsBs)) 6 

g *> 2a r(i-|) r2 V-l 1 2« 2 -l r(2-f)r(i-f) 



2(27r) d / 2 Vra 2 / 2 (2vr) d / 2 r (3 - d) 



(2tt)<V 2 I 2iU 2 



vn 1 ^ r 2 



ma 



' 2tt 



C + ln 



am 



V 2 y r(3-d) 

C = 0.577215... 



(d-2) 



Then, the functional £/ in (27) can be represented in the form 



ig(*B f ) . 



= yt/ + ^($)?7 1 + t/ J [$], 

where 

Q(g 2 a, A) = -g 2 N{a) - g 2 aA , R (g 2 a) = J da[£\ e~^ w ^ , 

U = [Da [l - e Q ^ a > A lR (g 2 a)] , U x = [Da a e Q ^ a ' A ^R (g 2 a) 



Ui[§\ = - [dx [ Da e Q ^ a ' A) [ da[£] e ~^ :W ^ : e~ 



-ig($B ( ) 
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The optimal value of parameter A should be obtained from the condition of elimination of 
the total linear terms over $: 



ig ($) 



D + A - -U x 



= 0, 



i.e. this equation reads 



A — —- I Da a 



The vacuum energy in the lowest approximation looks 

E = -g 2 (^-A 2 + D A - D^j + hj Q = -^A 2 + W , 
W Q = ^j Da{e Q ^ A) - 1 - Q(g 2 a, A) + e Q ^ a < A) [l - R(g 2 , a)]} , 
where divergences are completely removed as follows: 

-g 2 D A + g 2 D 2 - l - jDa Q(g 2 a, A) = . 
After complete removal of divergences we can put d = 2. Then, 



A 



oc 



da am' 2 



Ana 



1 - e Q{ - 9 



oo 



da am' 2 

e 2 



I Ana 2 
o 

Going to new re-scaled variables: 



'' a > A) R(g 2 , 

Q(g 2 a, A) _ x _ Q( g 2^ ^ + e Q(g 2 a, A) . ^ _ R ^ a )) 



with 



we obtain 



P = ® r , f (/?) = V® v( T ), k = -7=, 



N(t) =t[C + ln(t)] , a = — t 



g2a h + 
s = —— = at , 



h 



Z[g] = e- VE = e- VEo J D$ = e~ VE ° e~ VE ^ 

where the leading-order term for the vacuum energy is 



m 



and the remaining higher-order correction reads 

E corr = -±\nJ D$ e-^'W 
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(28) 



!1L- ) Eo = £ = -h { B 2 + [ - e-* 1 \e Q M - 1 - Q(ht) + e Q ^ (R(ht) -1)1 (29) 

07T / j J t ht L J 



(30) 



The shift parameter B is governed by the equation: 

1 r dt 



B 



t 



e Q(ht) R(ht) - 1 



(31) 



where 



Q(ht) = -ht[B + C + \n(t)} , R(ht) = J da[r]} e 



ht:Wh 



-±fdr vHt) 

do-[i]] = 5r] ■ e , 



da[r)} :W[rj\ : a = , 



± 

: W[rj\ : a = JJ dr x dr 2 



d( l . iq(v(ri)-v(T2)) . ^(in-ral-Cn-ra) 2 ) 

(2tt) 2 



The Green function for the kernel of da[rj] is 

\t — s\ t + s 



K(t,s) 



+ 



-ts 



Eqs. (p8|) , (p9| ), (|30|) and (|31~1) are basic in our consideration and they completely define the 
ground-state energy of the system at any given coupling. For general coupling these equations 
should be evaluated by numerical means, but we are able to solve them explicitly in the weak- 
and strong-coupling regimes. 



A. Explicit weak-coupling solutions 

i. First, we calculate the exact perturbation solution up to the order 0(g 6 ) as follows: 

e -VE( g)= l + 1 ^ 4 j y J d 2 xD A {x) + 0(5 6) 

So 4 V °r 3« 4 V 7 

= 1 + pfjw h « K ^ + °^ = 1 + jk^ 8^> + °^ • 

E -id^ 2 i a3)+0{3e) - h2 ^ 2 i <{3)+0<h3h 

where K Q (u) and ((x) are the Bessel and the Riemann function, correspondingly. 

ii. Now we evaluate E corr (g). We obtain the explicit weak-coupling solution to the higher- 
order total correction as follows (see Appendix C): 

n 4 7 m 2 7 

E co Ag) = -ttAt-, ~ C(3) + 0(g°) = ~h* — L C(3 ) + 0(h") . (32) 
{2Tx) 6 m z 8 57T 8 



iii. We have E = E +E corr by the very definition (28). Therefore, the explicit weak-coupling 
solution to the leading-order energy reads: 

n 4 7 m 2 7 

E = -7^r- 2 j Cm + 0(g e ) = -h 2 — - C(3) + 0(h*) . 
{liiym^ 4 Sti 4 
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Therefore, the following relations take place (within the given accuracy) 

E corr (g) = ~ E (g) = - E(g) . 

One can note that our leading-order energy E underestimates (i.e., 2/3 of) the exact energy 
E. The gap of one third is compensated by taking into account the higher-order energy term 
(|32|). As we will show below, our leading-order approximation becomes better as g grows and, 
approaches the exact solution for g — > oo. 



B. Exact strong-coupling solutions 

However, our method is designed to investigate the strong-coupling behaviour of the con- 
sidered system. Below we demonstrate explicit and analytic solutions of Eqs. fl28|), (|29| ) and 
(31) that leads to the exact strong-coupling answer. 

We note that for h — > oo the factor R(ht) gives correction ~ e o(\n(h)/h) w j^] 1 reS p ec t to 
the main exponentials in (|29"D and ([31]). Therefore, one can neglect it considering the main 
asymptotical values for B and E (for more details see Appendix B). 

For h — > oo the integrand in (|3T|) behaves a very sharply expressed Gaussian exponential 
centered at point t m . The exponent function is 

f(t) = -t- \n(t) - ht[B + C + \n(t)\ . 

Then, with the use of the "saddle-point" method we estimate 



B = \ J dtexv{f(t m ) - \ t 2 [-f"{t m )}] = \^U^)U^pu^y ( 33 ) 

where the maximum point t m obeys the conditions 

f'(t m ) = 0, f"(t m )<0. 

By resolving these equations we obtain (see Appendix D): 

B = \n(h) -C -2 + 0(1/ ln 4 (/i)). (34) 

Note, this asymptotics is reached very slowly, e.g., for h = 10 20 we obtain B = 0.93265 • \n(h) 
instead of B = ln(/t). Substituting (|5^ ) into ( ]2TS| ) and by neglecting terms vanishing as 0(1/ h) 
we obtain (Appendix D): 

S = -hi\n 2 (h) +3(C + 2)/iln(/i) +0(h). (35) 

Particularly, for h = 10 20 this asymptotics provides Sq = —2.82 • 10 23 while the numerical result 
is Sq = —2.81 • 10 23 . We see that this asymptotical behaviour becomes exact for very large 
coupling. 
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V. CONCLUSION 



We have represented a path-integral technique suitable to evaluate the quartic self interac- 
tion in the strong-coupling regime. The leading-order approximation is obtained in a relatively 
simple way and the remaining corrections can be estimated. A simple version of our method has 
been tested in the examples of a plain quartic-exponential function and the well-investigated 
anharmonic oscillator in Quantum Mechanics. We see that our technique effectively isolates the 
main contribution for arbitrary coupling. Then, we have applied this technique to the solution 
of the ground-state energy in the superrenormalizable scalar theory <f) A (x) in two dimensions. 
Hereby, the auxiliary appearing complex functional does not represent any difficulties for this 
technique. Our leading-order approximation becomes much accurate as the coupling g increases 
and for g — > oo it coincides with the exact asymptotics. We can conclude that our technique 
may effectively isolate the main contribution of the strong-coupling regime even in theories 
with divergencies and complex functionals. 
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TABLES 

Table 1. Approximate and exact numerical results for E(g) and 
G(g) in zero dimension. 



9 2 /2 


E (g) 


Ei(g) 


E{g) 


G (g) 


Gi(g) 


G(g) 


0.01 


0.02783 


0.02626 


0.02629 


0.90518 


0.90650 


0.90653 


0.1 


0.18027 


0.15267 


0.15361 


0.60786 


0.61378 


0.61553 


0.2 


0.27274 


0.22892 


0.23000 


0.49492 


0.50023 


0.50312 


0.5 


0.42431 


0.35925 


0.35993 


0.35854 


0.36211 


0.36596 


1.0 


0.55590 


0.47747 


0.47758 


0.27256 


0.27489 


0.27884 


2.0 


0.69826 


0.60939 


0.60890 


0.20316 


0.20461 


0.20823 


5.0 


0.89861 


0.79988 


0.79874 


0.13478 


0.13553 


0.13840 


10. 


1.05688 


0.95301 


0.95150 


0.09766 


0.09811 


0.10039 


100. 


1.60679 


1.49423 


1.49209 


0.03217 


0.03226 


0.03313 



Table 2. Comparison of results for E (g) , E osc (g) and 
Enum(g) in one dimension. 



9 2 
2 


E (g) 


Eosc(g) 


Enum(g) 


0.1 


0.06434 


0.05938 


0.05915 


0.5 


0.22666 


0.19697 


0.19618 


1.0 


0.35522 


0.30490 


0.30377 


10.0 


1.17291 


1.00778 


1.00497 


50.0 


2.30990 


2.00461 


1.99971 


100.0 


3.02802 


2.63759 


2.63138 
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Appendix A. 

The Green function in the presence of an external field reads: 

G m (t,1?\$) = —i ^ zn-rSCt-t?) 

v ' 1 ; D- 1 + m 2 + 2ig$(x) v ; 

s s 

°r n fdi-i-arh) ~ 2 W f dr *(t(r)) 

= J ds e -dW) . T J UmJ 9 { KK)) 5 (t-t>). 



Then, 

T e° o <S(* - f 7 ) 

s s s 

, - J*dr ^(r)+2 / ^(t) d _2i fl / dr $(t(r)) 

= Du e o o o _ 

, -/*" v 2 (T)-2ig f dr $( t+2 j dr' u(t') ] / ? 

= Du e o o V * JsU-t' + 2 drv[r) 

s s 

r °r Hs ~ sv l~S dr ? 2 (^)- 2 ^ / *(t+2(s-r)^o-2e(r)) 

= D£ du J- e ° o j(£ _ t' + 2sz/ ) 

— oo 

1 («-*') 2 /■ P{T)-2igfdT $(tl+t'(l-i)-2^(r)) 

e • / D£ e , 



«0)=£(s)=0 



where 



i/(r) + £(0) = £( S ) = 0. 



Appendix B. 

The following effective approximation takes place: 

R( s ) = J da[rj\ e~ sW = ( e- sW ) « e" sW • cosh (V^ 2 ) - (W^) 2 ) , W =: W[r}\ : a . 

Due to normal-ordered form of FT, one finds (W) = 0. Then, 

R(s) ss cosh(sw) , u? 2 = y dcrfa] (: Wfo] : CT ) 2 . 

For /i — > oo its correction to the shift parameter B = ln(h) is ~ 0(1) and to the leading-order 
energy E = — (3/2)Mn 2 (/i) - ~ 0(h). Therefore, its influence is negligible and we just drop 
the factor R(ht) by considering the large h asymptotics. 

Appendix C. 

By expanding the interaction functional 

U l m=g 2 U 2 l<5>]+g 4 U 4 m + 0(g 6 ) 
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and taking into account its normal- ordered form, we obtain 



E C orr(g) = - J^ JjDaDP JJda^daW JJdxdy Jd<$> : {^f xa :: (*B v ) 2 yfl : + 0(g 6 ) 



The functional averaging over field <£> results in 



Jd<$> : :: {QB^ := 2 ^cftdu ffdsdw J-^- j-^- ^(fc+rt^-^+^KW-K^l+^KH-^-)] 



Further we use the following plain relations: 

yka; e l{k+p){x - y) = V • (2tt) 2 • 5 2 (k + p) , 
Jda[i] 

a 

JJdt ds e -£*.(l*-l> = 2a Jdt e-^-W (l - £) , 



e <fcKW-«-)l = eX p 







fjpk_ -!^.[F a (t)+Fp(8)] _ J_ fj p -§[F a (t)+F fl ( S )] 



Then, 



(2tt) 2 





^ 1 00 00 



S corr (<?)= -grf^i JJdxdy(l-x)(l-y) jdu JJdad/3eS ^ 2+ux (i-*)]-f K+^Ml 
^ ' 00 

# 4 fh j r-t \/i Jn[x(l -x)] -ln[y(l - y)] 



fcdy(l - x)(l - y ) ^- J 



lQir 3 m 2 JJ x(l — x) — y(l — y) 

91 I as). 



(2Txfm 2 8 

Appendix D. 

The maximum point t m of the sharp exponential function 

f(t) = -t- \n(t) - ht[B + C + \n(t)} 

is dictated by the conditions 

f(t m ) = -1 - — - h[B + C + 1 + ln(t m )] = , f"(t m ) = < . 

This results in 



B = \n(h) - ln(0 - C- l- i- i = , e = ht m > 1 . 
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Hereby, 

The solution is 



lim = +1- 

n^oo 



B = \n(h) - C - 2 + 0(1/ ln 4 (/i)) . 
For the energy, we rewrite it in a more convenient form 



dt 



where 



S = -hlB 2 -Jj e"* - $(M)] , 



«A „ . ,.s „ -\n+l 6 



F(t) = /^e- s + ln(t) + C = E(-ir 

/ n=l 



n\ n 



and 

$(^) = _L | e -w[B+c-in(h)+in(«)] _ i + ht[B + C - ln(/i) + ln(M)]} 



ht 

Therefore, 

oo , 



and 

eft 



[ 1 e -t $(^) = I + 0(ln(/i)) . 



Substituting these solutions to the energy, one obtains 

3 



S = --h\n 2 (h) + 3(C + 2)h\n{h) + 0(h) . 
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